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1. Introduction 

According to Shannon's information theory [I], the only rigorous measure of the 
uncertainty or lack of information associated to a continuous random variable X with 
density function p{x), X e R^, is the entropy 

H{X) = - J p{x) log p{x) dx . (1) 

In particular, when p{x) is the single-particle probability density for position of 
a quantum system, H{X) is the only rigorous measure of the uncertainty in the 
localization of the particle in position space. The momentum entropy H{P) can be 
defined likewise from the single-particle density of momentum 7(p). In the simplest 
case of a single-particle system described in position space by the wave function ip{x), 
we have that p{x) = and 7(p) = where the wave function in momentum 

space is the Fourier transform ofip{x). The sharp inequality [2l [3] 

H{X) + H{P)>D{l + \og7i) (2) 

places a nontrivial lower bound on the sum of the uncertainties in position and 
momentum, so it provides a quantitative formulation of the position-momentum 
uncertainty principle. Using the variational inequality that relates information entropy 
and standard deviation for an arbitrary D-dimensional random variable [HIS], 

ff,^)<«fl + ,ogH^V (3) 



2 V D 

the entropic uncertainty relation (I2l) leads to the well-known Heisenberg uncertainty 
relation 

AXAP > ^ , (4) 

which proves the former to be stronger than the latter. 

For many important quantum systems, such as D-dimensional harmonic oscillator 
and hydrogen atom, the calculation of position and momentum information entropies 
involves the evaluation of integrals of the form 

E{Pn) = - {pn{x))'^ log {Pn{x) f Uj{x)dx , (5) 
J a 

where {pn{x)} denotes a polynomial sequence (degp„(x) = n) orthogonal on [a,b] C M 
with respect to the weight function uj{x). During the last decade there has been an 
intense activity in the study of these integrals, motivated not only by their relevance to 
quantum physics but also by their close relationship to other interesting mathematical 
objects, such as the L^-norms or the logarithmic potentials of the polynomials Pn{x). A 
survey on the state-of-the-art in this field up to year 2001 can be found in [4]. 

The calculation of the entropic integrals E{pn) is generally a very difficult task, 
and in most cases only asymptotic results for large values of n are known [4]. In fact, 
since all the zeros of p„ are simple and belong to (a, 6), when n is not very small even a 
numerical computation of E{pn) poses serious difficulties due to the strongly oscillatory 
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behaviour of the integrand in ([5]). In this respect it is worth mentioning Ref. [5], which 
presents an efficient algorithm for the numerical evaluation of E{pn) in the case when 
the interval (a, b) is finite. 

Closed analytical formulas for E{pn) are only known for a few particular cases of 
the Gegenbauer or ultraspherical polynomials We recall that these polynomials 

are defined as (see, e.g., |6l Sec. 4.7]) 

C^^^Hx) = 7^Pi"^'"^\x) , (6) 

where {a)n = r(a + n)/T{a) denotes the Pochhammer symbol and P^'^\x) are Jacobi 
polynomials, 



For A > — |, Gegenbauer polynomials form an orthogonal sequence on the interval 
[— 1, 1] with respect to the weight function wx{x) = (1 — x^)'^~2 ^ 

; c;?>(.)e(.)(i - .Y-i d. ^ QI^^„,.„ . (8) 

The information entropies of Gegenbauer polynomials, on which we focus in the present 
paper, are thus given by 

i?(CW) = - j'^ iCi'\x)f log iCi'\x))\l - x'f"'^ dx . (9) 

The integrals E{Cn^^) are especially relevant in the case when A is a non-negative 
integer or half-integer number, due to the relationship between the corresponding 
Gegenbauer polynomials and (hyper)spherical harmonics. As a consequence, these 
integrals appear in the calculation of the angular component of information entropies 
in both position and momentum space for any D-dimensional (D > 2) quantum- 
mechanical system with a central potential, such as the isotropic harmonic oscillator 
or the hydrogen atom (radially symmetric Coulomb potential) [H [3, El [9]- They also 
control the radial component of the information entropy in momentum space for the 
D-dimensional hydrogen atom [11 [3, [8]. 

Instead of using the standard definition of Gegenbauer polynomials, it is often more 
convenient to work with the polynomials 

e?'w = (^^^?c«^'w. (10) 



A(2A) 

which are orthonormal on [—1, 1] with respect to the probability density 

^^A(x) = -g^±^(l-xy-K (11) 

The corresponding entropies, 

E(aW) = - j\ci'\x)]Hog[Ci'\x)rw,{x)dx , (12) 
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are related to E{Cn^^) by the formula 

which readily follows from the previous definitions by taking into account the 
orthogonality relation (JHj). 

The simplest particular cases of Gegenbauer polynomials are the Chebyshev 
polynomials of the first and second kind, 

For both of these families, information entropies can be computed in closed analytical 
form, the results being [71 [8] 

f if n = , 

I log 2 — 1 II n 7^ , 

E(U.) = . (16) 



In the A = 2 case, it was first proved in [TO] that 

F(r^'^)- 1n.r ^(^ + ^^ ^ n{n^ + 2n-l) 2 T^^IQ 

^\ n + 3 J (n + l)(n + 2)(n + 3) ^(^ + l)3(n + 3)3 T^'+^IO ' ^ ^ 



+ l)(n + 2)(n + 3) ^(^ + i)3(^ + 3)3 T^'^^IO 

, ^ + ^ , (18) 
v/(^ + l)(n + 3)' 



where 



and this result was later simplified to [TT] 

(19) 



£;(c(2)) = -iog( 



3(n + l)\ - 5n^ - 29n - 27 1 /n + 3^"+^ 



n + 3 J {n + l){n + 2){n + 3) n + 2\n + l 

In the same work [TI], it was also obtained the following generalization of (fTTl l to 
arbitrary integer values of the parameter, A = / G N: 

where Sni and r^; are known constants depending only on n and /, the auxiliary 
polynomials P and H are defined from the sequence {Pk} (degP^ = k) generated by 
the recurrence relation 

Pi,^^{x) = {21 -2k- 3)xPk{x) - {n + k + l){n + 21 - k - - x'^)Pk-i{x) (21) 

from the initial values P-i{x) = 0,Po{x) = 1 through the formulas 

21-2 

P{X) = P2l-2{X) , H{X) = J2i-'^yPs-l{x)P2M{^) , (22) 

s=Q 
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and {j = 1,2, ... ,21 — 2) denote the zeros of P. The explicit expression of the 
polynomial P was later found to be 



in + l)n\ j^^{l-n- 1)^(1 + n + l)^fi\^ ^ ' 

Regretfully, (l20l) is not easy to use in practice. Furthermore, it is not a completely 
analytical formula save for small values of / since, as we readily see from (1231) . the zeros 
of P have to be determined numerically when / > (j^. 

As first pointed out in |T3], the entropy of Chebyshev polynomials of the first and 
second kind can be easily computed by direct calculation of the corresponding integrals 
by using the well-known trigonometric representations 

T„(cos^) = cos 720 , UnicosO) = ^^^(^ + ^)^ ^ (24) 

with X = cos 6. Motivated by this observation, in the present paper we aim at evaluating 
the entropic integral E{Cn^^) for general values of the parameter A using representations 
of the same kind for the Gegenbauer polynomials. 

We begin by collecting, in Section O the trigonometric representations of 
Gegenbauer polynomials that will be used later on. Our approach is developed in 
Section [3l where we show that it enables us to find completely analytical expressions 
for E{Cn^^), in terms of finite sums, whenever A G N. The new results obtained for the 
information entropy of Gegenbauer polynomials of integer parameter are summarized 
in Section [H Finally, in Section [5] some concluding remarks are given and several open 
problems are pointed out. 

2. Trigonometric representations for Gegenbauer polynomials 

The most widely known trigonometric representation of the Gegenbauer polynomials is 
(see e.g. p] p. 302]) 

n n 

(cos e) = Y^ ^V^"-''"^' = E - 2"^)^ ' (25) 

m=0 771=0 

where 

MX) 

" m!(n-m)! ' ^ ^ 

Another representation, due to Szego [6l[l5], is 

(A) oo 

CW(cos0) = ^-^^2A^E"Ssin(^ + 2z/ + l)0, A>0,A^N, (27) 

where 



(,) ^ 2^-^^r(n + 2A) (,) ^ (1-A).(n+1). 

" r(A)r(n + A + l) ' "'"^ i/!(n + A + l), ■ ^ ^ 

I Likewise, the general expression of E{C^^) given in [9] is not completely analytical save for small 
values of n, since it is expressed in terms of the zeros 
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At first sight, this representation seems to be less useful than the previous one, because 
it contains infinitely many terms. Moreover, it is supposed not to hold when A G N. 
However, it is not difficult to prove that the validity of ( l27ll extends to the case when A 
is a positive integer. 



Proposition 1 The Szego representation [27]) holds true when A G N. In this case, it 
reads 

(A) A-1 

Cl^) fees 6) = + + 1)^ . (29) 



sm I 



n2A- 



u=0 



Proof If A G N then = when z/ > A, so that (ET]) reduces to ((29]). We will prove 
this equality by induction on A. When A = 1 (l29ll is obviously true since it reduces to the 
second equation in f l24l ). the well-known trigonometric representation for the Chebyshev 
polynomials of the second kind. Now, assume that the result holds for A = m — 1 
(m G N). We take advantage of the following recurrence relation for the Gegenbauer 
polynomials [H Eq. (4.7.27)], 

nC^'\x) = {2X + n- l)xC^^\{x) - 2A(1 - x^)^^!^^^), (30) 

which in trigonometric form {x = cos 6) can be restated as 



CW(cos^) 



(2A + n - 1) cos^Ci+7^^(cos^) - {n + 2)0^^+2'' (cos 6) .(31) 



2(A- l)sin2^ 

Using this formula for A = m and substituting (i29l) on the right-hand-side we arrive at 



2ct'> 



sm( 



\2m-l 



m— 2 



cos 9 ai^+i^ sm{n + 2u + 2)9 



u=0 



n 



m~2 



n + m + 

(m) 



I • n\2m—l 

{sm9) 



Y sin(n + 2u + 3)9 



u=0 
'm-2 



m-2 



Y Oi'Z'+l sm(n + 2u + i)9 + Y "fc'i^ sin(n + 2u + 1)9 



2{n + 2) 



,u=0 
m-2 



u=0 



I sin(n + 2z/ + 3)9 



n + m + 



u=0 
m—1 



sm( 



\ 2m- 



a 



(m-l) 



n+1 "I" ^i/. 



(m-l) 



n+1 



u=0 

X sin(n + 2z/ + 1)6' 



2{n + 2) 
n + m + 1 



-a 



(m-l) 
1^-1, n+2 



(32) 

where in the last step we have used that = and = whenever u < 0. 

A straightforward calculation shows that 



(m— 1) (m—1) 



2{n + 2) 
n + m + 1 



-a 



(m-l) 
u-l,n+2 



a 



(m) 



(33) 



and (l29l) is thus proved to hold also for A = m. □ 
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The fact that the sum in ( !27ll terminates after a finite number of terms when A G N 
suggests that Szego's representation may be useful to evaluate the entropy of Gegenbauer 
polynomials of integer parameter. Accordingly, in what follows we shall assume that 
A G N unless otherwise indicated. 

3. Evaluation of the entropic integral 

With the change of variable x = cos 9, the integral ([91) takes the form 

E{Ci^'>) = - r (cos 6) f log (CW (cos 6) f sin^^ 6 dO . (34) 
Jo 

Using Szego's representation (l29l) for one of the two Gegenbauer polynomials in 
{Cn (cos 6')) , (1341) can be rewritten as 

E{Ci^') = g - Jii\,.) , (35) 

where 

jW := / C(^)(cos^) cos(n + 2u)e\og {Ci^\cose)fde . (36) 
Jo 

Now, using the standard representation fl25l) we have that 

^^(cos^) cos(n + 2u)e = -J2 dm,n cos 2{m + J^W + -J2 ^m,lcos2(n - m + z/)^ .(37) 

m=0 m=0 

Taking into account the symmetry property (im,n = d^n-m,m which readily follows from 
the explicit expression of the coefficients (im,n, the previous equation simplifies to 

n 

CW(cos^) cos(ri + 2u)e = ^ c/^^)„cos2(m + u)9 , (38) 

m=0 

SO that 

4J] = E<1 rcos2(m + z/)^log(CW(cos^))'d^ . (39) 

771=0 

Defining the integrals 

4'i:= rcos(2me)log(CW(cos^^))'de, (40) 

^0 

from ( l35l l and ( l39ll we find that E{Cn^^) is given by 

EfC^^)) = --c(^) Va(^) Vrf(^) f/^i^ -J^^l ^ (41) 

u=0 m=0 

An alternative expression for the entropic integral E(C-n^) which turns out to be 
more convenient in practice can be obtained by noticing that 

^ A— 1 ( n n+1 

EiC^^^A = -- c(^) V a(^) I V c/(^) J^^^ - V I^^^ 

J 2 / J u,n I / J "'m,n u+m,n / j "'m—l,n u+rn 

u=0 l.m=0 m=l 

X-1 ( n ^ 

- _ic(^) Va^^N V ('rf(^) -rf^^^ -rfW/^) I (42) 

1^=0 lm=l J 
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According to f l26ll . dm^n = when — A<m<Oorn<m<n + A. Let us restrict 
initially to the case when A 7^ 1, so that = and d^n+in — 0- This allows us to 

write the previous formula in the more compact form 

A-l n+1 



7-(A) 



u=0 m=0 
A-l n+l+u 

= _ic(^) V V a(^) ("rf^^^ -rf^^^ ^/(^) (4-]) 

2 n / y / ^u,n y^m—u,n "'m—u—l,n J m,n • V ^/ 

i/=0 m=v 

Usina; = when —A < m < as well as when n < m < n + X, we 

can extend the lower and upper limits in the inner summation to 1 and n + A — 1, 
respectively, provided that the terms m = and m = n + X are treated separately. Thus 
we find that 

-, / n+A-l \ 

^^r'(A)A__icW r^{A)^(A)j{A)_ (A) ^(A)j(A) , ^(A) j(A) /^^n 

^K'^n J — 2 I 0,n"0,n-'0,n "A-l,n"n,n-'n+A,n ' /-^ Pm,n''m,n I ' \^^) 

\ m=l / 

where 

A-l 

ym,n / J ^u,n \^m—u,n "'ra—v—l,nj • V "/ 

It can be seen that (l44l) also holds when A = 1 by noting that in this case its right-hand 
side coincides with that of fl42l) . 

In order to apply (l44l) . we need to evaluate the integrals /m,n with < m < n + A. 
This goal can be achieved by means of complex integration techniques, which enable us 
to obtain the following result. 

Theorem 1 For A G 

/™ = 2.1og(<^) (46) 



and, when m > 1, 

(47) 



m ^ (2m)! rf^^™ 1 iogZ^«.,nl^ J 



i/=0 



2 = 



Proof Taking into account that Cn^\—x) = {—iy"'Cn^\x), (i40l) can be written as 
= 1 cos(2m^) log |CW (cos 9)fd9 

cos{2m9)\og\C^^\cos9)\d9 

/ cos(2me) log |e^"^C(^) (cos 0) | d^ , (48) 
Jo 



2lT 
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where in the last step the factor e'"^^ has been introduced for later convenience. Using 
the Szego representation (l29l) for the Gegenbauer polynomial Cn'^^ (cos 6') , the previous 
equation reads 

f27r 



JW„= / cos(2m^)log 
Jo 



^^5^aWsin(n + 2z. + i: 
sm 6 

i/=0 



de . (49) 



We will compute the integral 

'"27r / ^(A) i„6l ^-1 



cos(2m^) log I -^^Q E sin(r^ + + 1)9 \ dO , (50) 



u=0 

r(A) 



whose real part equals Im 'n- Introducing the change of variable z = exp(i6') we arrive at 

J\z\=l ^ 



where 



= -"c^f ^ (^^) = ^i'^2^''^(-i)' ^^^°""";;_,2)2A-r — ' ■ (52) 

The singularities of the integrand are z = 0, which is a pole of order 2m + 1, and all 
the zeros of q{z), which are branch points. If {xnj}^=i denote the zeros of 
which are known to be simple, real and located in (—1,1), then the zeros {zn,j}^Zi of 
the function q{z) are 



z„j+Hzp| = exp(i arccosxnj) = Xnj ± i-^1 - x^j , j = 1, 2, . . . , n . (53) 

This means that the {znj}'j=i are all located on the unit circle, which can also be seen 
from the fact that z = exp(i arccosx) maps (—1, 1) onto the unit circle. Therefore, the 
integrand of (ISTl) has 2n branch points located on the contour of integration. To avoid 
this difficulty we consider the same integral along the slightly different contour F (see 
Figure [T]), which is also closed. Notice that the logarithmic branches can be chosen to 
go from the branch points to the exterior of the unit disk, so that F does not cross them. 
Since the only singularity inside F is 2; = we now have 

r _|_ ]^ / ^4m _|_ ]^ \ 

^2m+i logg(^)d2 = 27riRes(^ ^^^^^ logqiz), z = j . (54) 
The integral along F can be decomposed as 

(55) 



7=1 \ "^J "ii 



where ej denotes the arc of circumference of radius e that surrounds the branch point 
Znj and 7j denotes the arc on the unit circle that connects ej and e^+i (72n connects e2n 
and ei). Parameterizing z = Znj + we find that 

/•;z^™ + l, f (Zni + e&^Y'^ + l 



/(^ . _|_ ppio^m _|_ 1 
J)2mti ^(--.- + ^^^') '^^''dO , (56) 
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Figure 1. Integration contour used to avoid the branch points on the unit circle. 



where we have used that xlogx — >^ as x — 0. Thus, taking the limit e — > in (l55l ) we 
conclude that 



\ogq{z)dz 



(57) 



Taking into account that the residue of a meromorphic function h{z) in a pole zq of 
order 2m + 1 is given by 

1 d^m 



Res{h{z), z = Zq) 



(2m)! ciz^m 
use of (iniD and ^ into ^ leads to 



((Z - Zo)''"+'/^(^)) 



2 = Z0 



j(A) 



[(/"^ + l)logg(z 



2=0 



(2m)! rf^^m 

In the case m = 0, the previous equation reduces to 

Xii; = 2vrlogg(0) = 2vrlog [cl^'h^'-^-l) 

so that 



"A-l,n / ' 



'0,n 



27rlog 



^{A)22A-2^(A) 



277 log 



(A)n 



(58) 



(59) 



(60) 



(61) 



which proves the first part of the theoren|§|. On the other hand, if m > 1 then we 
readily see from (i59ll that Xm,n € M, so /m,n = Xm,n- Furthermore, in this case the factor 
(z^™ _l_ ^j^g right-hand side of fl59ll can be omitted, since at z = its value equals 
unity while all its derivatives do vanish. We thus find that 

^A-l ^,(A)^^2?i+2A+2i^ _ ^2A-2i/-2^| 



r(A) ^ ^ ^ 

"^•^ (2m)!rfz2m 



log 



1-^ 



2^2A-l 



(62) 



2=0 



This part can also be proved using the mean value theorem for harmonic functions (cf. |7i Sec. VI]). 
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and f l47ll follows on noting that 
(log(l - z')) 



2m 



dz 



2m 



E 



k=l 



(2m)! 



m 



(63) 



when m 7^ 0. □ 

The fact that the Szego representation ( l29l) has a finite number of terms plays 
an essential role in the proof of Theorem [TJ Although we are mainly interested in 
evaluating the integrals /m,n when A G N, it is worth pointing out that these integrals 
can be calculated in a similar way for all possible values of A, provided that we use the 
standard trigonometric representation (l25ll instead of the Szego representation for the 
Gegenbauer polynomial inside the logarithm. This generalization is contained in the 
next theorem. 



Theorem 2 For XeR, \> 



1 

' 2> 



-'O.n 

and, when m > 1, 



2lT log 



(A), 



j(A) 



TT 



d' 



'2m 



(2m)! d^^m 



(64) 



(65) 



2=0 



Proof We proceed as in the proof of Theorem [H but now we use the complex form of the 



standard trigonometric representation (l25ll for the Gegenbauer polynomial Cn^^ (cos 



m 



Thus we arrive at 



j(A) 



2m 



m,n 



^4m 



0=0 



(66) 



z=0 



(2m)! 

from which (l64ll and (l65l ) readily follow. □ 

In order to carry out the sums in (i44l ). the next step is to obtain closed formulas 
for the derivatives Im]n with 1 < m < n + A. Despite its greater generality, Theorem 
[2] turns out to be less useful than Theorem [1], because ( l47l ) expresses the integrals in 
terms of the logarithm of a polynomial that has 2A terms, while in (l65l) they are given in 
terms of the logarithm of a polynomial with n + 1 terms. As we shall see, the difficulty 
in obtaining a closed formula for the derivatives of such functions increases with the 
number of terms in the polynomial. Therefore, if we want an expression of /m,n for a 
fixed value of A and any n G N Theorem [U is more helpful, particularly for small values 
of A. 

In the case A = 1 we readily notice from f l47l ) that, ifl<m<?2 + l, then 



J(i) = - + - 
'"■^ m (2m 




2 = 



2 = 



(67) 
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When A > 2, the polynomial inside the logarithm has more terms and the above trick 
does not work. However, we can obtain closed formulas for the derivatives in (l47l) by 
means of Faa di Bruno's formula for the derivatives of the composition of two functions, 
which states that (see e.g. [T6] ) 



m 



A;=0 fci,fe2,...,fc„ j=l 



(68) 



where the inner summation is extended over all partitions satisfying 

h + k2 + \- km = k, ki + 2k2-\ h mkm = m . (69) 

This formula enables us to find explicit expressions for Jm.n with A > 2, as stated in the 
following two propositions. 

Proposition 2 In the case A = 2, when 1 < m < n + 2 



/(2) = L 



n + 3 
1 



n 



n + 3 

+ — Om n+2 

n + 1 



(70) 



Proof In this case, application of Faa di Bruno's formula (l68l) to the derivatives in (l47l 
give 



2m 



dz 



2m 




2m 



d^ 



2 = 



fc = l 



X E n 

On the one hand, for /c > 1, 



2m f 



al^(z2-+4+2z^ - ^2-2.)) 



z=0. 



{3\pk,\ 



so that 



d^ 
dz^ 

d^ 

dz^ 



(log z) 
(log z) 



{-lY+\k-l)\ 



(71) 



(72) 



(73) 



On the other hand, all the derivatives of the polynomial in (uTl) vanish at z = except 
when j = 2 and j = 2n + 4, so we must set kj = if j 7^ 2 and j ^ 2n + 4. Conditions 
flHO]) then read 

^2 + k2n+i = k, 2k2 + (2n + 4)A;2„+4 = 2m . (74) 

Since A;2 and A;2„+4 are non-negative integers, these equations only admit the solution 
fc2n+4 = 0, ^2 = = m when m < + 1, while in the case m = n + 2 we have to add 
the solution fc2n+4 = A; = 1, A;2 = to the previous one. Therefore, fFfTl) simplifies to 



d 



2m 



dz'^'^ 




(2) / 2n+4+2i/ _ ^2-21/ \ 



-(2m)! 



2 = 



a, 



(2) 
0,n 



(2) 



+ 



a 



a 



(2) 

0,n (- 

'-'m,n-|-2 



1,71 



,(75) 



and the result follows using the second equation in (1281) . □ 

II Notice that attending to ((69|) /c = corresponds to m = 0, so we can start the sum in fc from 1. 
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Proposition 3 For any A G N, A > 3, when 1 < m < n + \ 



a 



TC- 



(A) 
0,ra 



m 



a 



(A) 
A-l,n 



6^ 



m,n+X 



a 



TT 



(A) 
A— 3,n 



a 



m [m/3] [m/4] [m/(A-l)] 



(A) X 
A-2,nJ 



(A) 

A-2,n/ A;=l ^6=0 fe8=0 



'«2{A-l)=0 



(A) (A) 
'-'^A-l,n'-'^A-3,n> 



{k-l)\ 



{2k -m + Y.t=i rk2r+4y- {m - k - Y.'lJ^i^ + 1)^4)' 
A-l . / (A) ' (^'^ -^''2.' 



^A-3/ 



/ (A) \i-2\"^^ 
«A-l-i,n(«A-2,n)' 



i=3 



(A) 
*^A-3,n/' 



(76) 



where in the upper limits of the summations over k^, ks, ■ ■ ■ , k2{\-i) the square brackets 
denote integer part of the expression within. 

In particular, in the case A = 3, when 1 < m < n + 3 

(n + 4)(n + 5) 



where 



/(3) 



f{n) 



TC 



m 



5 - 23fJ(/H 



— TT- 



[n 



(n + l)(n + 2) 
+ l)(n + 5) + iv/3(n + l)(n + 5) 



5. 



m,n+3 ; 



(77) 



(78) 



(n + l)(n + 2) 

Proof In the general case (A G N, A > 3), application of Faa di Bruno's formula ([68 
to the derivatives in (l47l) and use of (rf2l) lead to 



2m 



A-l 



2m 



{k-i)\ 



2n+2A+2f _ ^2A-2!^-2\ 



u=0 
2m 



2m 



(2m)! 5^ 



fc=l ['^\^l,nJ A:i,fc2,...,A;2m i=l 



(A) 



E n 



z=0 

^2n+2A+2i/^ 



2 = 0J 



(j!)'=.A;,! 



2m 



/ ,,>r^ (k-iy. 



2m 



fc=l ['^\-l,nJ ki,k2,...,k2m i=l 

In the first term of the right-hand side all derivatives vanish at 2; = except when 
j = 2n + 2A, so that kj = whenever j ^ 2n + 2A and conditions (l69|l simplify to 

fc2n+2A = k, (2n + 2A)A;2n+2A = 2m , (80) 

which only admit the solution k2n+2\ = k = 1 when m = n + A. In the second term the 
derivatives that do not vanish are those with j even, 2 < j < 2A — 2, so that conditions 
(l69l) now read 

A-l A-l 

'^k2r = k, '^sk2s = m. (81) 



E n 



_ v^A-1 (A) 2A-2!/-2 



2=0. 



(j!)^.fc,! 



(79) 



r=l 



s=l 
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Equation (i79l) thus reduces to 



^2m / A-1 



dz 



2m 



u=0 



z=0 "A-l,n 



-(2-)'L-(Ar^ n — ik^i — • (82) 

Finally, we can further simplify the previous expression to obtain (l76ll by using 
conditions (ISB to write ^2 and /C4 in terms of the remaining indices, i.e. 

k2 = h + 2k8 + --- + i\- 3)A;2(A-i) + 2k - m , 

ki = -2ke -Sks (A - 2)k2(x-i) +m-k. (83) 

Notice that in (TfGll conditions (ISTll are guaranteed to hold because for the values of the 
indices that do not fulfill them we get the inverse of the factorial of a negative integer, 
which can be considered to be zero. We have changed the upper limit in the sum over 
k from 2m to m because when m + 1 < k < 2m conditions (jSTll are not fulfilled. 
In the case A = 3, (1761) reduces to 



2,n 



SO we need to evaluate a sum of the form 

v ,,, ,,, x^=E ^?^';■V^ '"^^^=^"E^""0^^ (85) 

{2k - m)\{m - k)\ JT^ jK"^ " 2j)! V J^^-J 

Using the summation formula [171 Eq. (5.75)] 

|(70^-(^)'"^(^)'" 

and the second equation in fl28l) . we find that 

- ffr^^-^ prr,^™i ^ (^ + 4)(^ + 5) 

^-''^ - m " m f^(^) + / (^) ] - ^(^ + i)(^ + 2) ^"'-^^ ' (^^^ 
which is equivalent to f TTTl) . □ 

4. Results for the information entropy 



Equations (1461 ). (i67l ). (1701 ) and (iTTj ) enable us to derive closed analytical formulas for 
E{Cli^) when A = 1,2,3. For A = 1, after substitution of the corresponding values of 
the constants (im,n, ci^'' and (see (1261) and (l28l)). ( l44l) reduces to 

^(Ci'^) = -^(/S-4^i,n), (88) 
which using ( l46l ) and ( l67l) immediately leads to 

=£(£/„) = ^(^-1). (89) 
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E{Ci'^) = ( + l)(n + 3)45 - + l)'^S2,n - 4 J2 ) ' (90) 

\ m=l / 

SO using ( l46ll and (FTOl l together with the well-known formula for the sum of a geometric 
series, 

x(l — x") 



E 



m=l 



1 — X 



(91) 



we find that 

) = - 



TT / , , , , , , , ra^ - - 29n - 27 
2(r2 +l){n + 3) log(n + 1) + 



n + 2 



(n + 3) 



ra+3 



^(n + 2)(n + l)"+i j ■ ^^^^ 
Recalling (fT3l) . (l89l) and (l92l) are readily shown to be equivalent to (fTBl) and (fT9l) . 
respectively. 

In the case A = 3, fl44ll can be rewritten as 
1 



E(C7f) 



128 



(n + l)(n + 2)(n + 4)(n + 5)4 

n+2 

n+l)2(n + 2)2/g3,„-12^ 



(3) 



mm 



6n + 7-2m2)/W^ 



m=l 



(93) 



Substituting (l46l ) and f FTTl ) into the above expression, we encounter again the geometric 
sum, as well as a sum of the form ^^m^x'". Using (l9Tll and the summation formula 
dH Eq. (5.14.9)] 



m=l 



r, ^ x(l + x) - x'^+^\(n + ly - (2n^ + 2n- l)x + n^x^] , , 
m X = ^ 77 , (94) 



(\-xf 



after a tedious but straightforward calculation we arrive at the following closed analytical 
formula for E{Cn ), which is a new result: 



E(C(^)) 



TT 



128 



2(n + l)(n + 2)(n + 4)(n + 5) log 



[n + l)(n + 2) 



rc" - 16n^ - 269^3 - 1200^^ - 2102n - 1250 



+ 



(n + 2)(r2 + 3) 



n + 3 

{n + l)(n + 5) + i^3(n + l)(n + 5) 



n+l 



(n + l)(n + 2) 



X 2n^ + 13n + 14 - i(ra + l)(n + 6) 



> + l)(n + 5) 



(95) 



When A > 4, combination of (l44l) and (TfGl) provide an expression for the entropy 
-E(Ci'^-') in terms of finite sums. For the sake of brevity, in (l44l) it is convenient to absorb 
the term corresponding to /^^\ „ into the sum over m by setting 



(96) 
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instead of using ( l45ll with m = n + X, which would give for I3^n+\n the value 

-ot^\li^n'^^n}n + «o!'i'^i+A,n- ^c thus havc that 

E{Cl>^') = - f A^id^log + + (2A - 1) X: 



n! / ' ' — ' m 

m=l 

n+A m [m/3] [m/4] [m/(A-l)] / (A) \™ / . (A) ^2 

W V V ■ ■ ■ V ^-3'" 1 _^^^A-2^ 

I (A) II (A) (A) 

m=l fe=l A;6=0 fc8=0 fe2{A-i)=o \'^\-2,nJ \ '^A-l,n'^A-3,n , 



(A:-l)! 



{2k -m + Y.t=i rk2r+4V- {m - k - Y!1=i{s + l)^2s+4 



X 

j=3 V V"A-3 



Unlike (l20l ). (l97l) is completely analytical for all A G N, which makes it suitable for 
symbolic computation. For instance, a Maple implementation of the formula enabled us 
to obtain the closed analytical expressions for E{Cn^) and E{Cn^), with 1 < n < 15, 
that are displayed in Tables [H and El respectively. In these tables we also provide 
numerical values of the entropies obtained from the exact ones, in order that the 
interested reader can compare them with those given by numerical algorithms such 
as that in 0. 



5. Summary and conclusions 

The problem of obtaining closed analytical formulas for the entropy of orthogonal 
polynomials is known to be very difficult, as displayed by the fact that in previous work 
on the subject formulas of this kind were only found for the Gegenbauer polynomials of 
parameter A = 0,1,2. Here we have presented a new approach to the calculation of the 
entropy of Gegenbauer polynomials, based on the use of trigonometric representations 
for these polynomials, which has allowed us to explicitly evaluate the entropic integrals 
by means of complex analysis techniques. Using this method we have been able to 
derive in a unified way closed formulas of E{C^^) for A = 1, 2, 3, the last one being new. 
Furthermore, when A > 4, A G N, we have obtained completely analytical expressions 
of the entropy in terms of finite sums, which easily provide exact values for the entropy 
using symbolic computation. The growing complexity in the formulas of E{C^^) as A 
increases serves as a clear illustration of the difficulties posed by the calculation of the 
entropy of orthogonal polynomials. 

When the parameter A is not a positive integer, the Szego representation ( l27l) of 
the Gegenbauer polynomial has infinitely many terms, so the same happens for the 
expressions (l35ll and (1411) of the entropic integral E{Cn^^). It remains open the problem 
of studying the convergence behaviour of these series, as well as that of summing up them 
analytically. It would be of particular interest to obtain exact analytical expressions 
for the entropy of Gegenbauer polynomials of half-integer parameter since, as already 
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Table 1. Exact and numerical values of the entropy E[Cn ) for 1 < n < 15. 



n 


Exact value 


Numerical value 


1 


119 

-77rlog(2) + ---7-7r 
' 240 


-13.685 


2 


105 , , , 580771 

;^7rlog 10) + ■— — —ttt 

8 ' 300000 


-88.862 


3 


75 , 95 
-— Trlog 20 + — TT 
2 ^ 21 


-338.714 


4 


5775 , ,^ , 4883222845 
64 ' 632481024 


-983.613 


5 


385 , 17355685 
2 ^ ' 1806336 


-2404.173 


6 


3003 , , , 6449434961 
8 ^ 1058158080 


-5206.005 


7 


1365 , 1396715852287 

77- TT log 120) , TJ" 

2 6^ ^ 139218750000 


-10296.556 


8 


75075 , , 24757176334716125 

TT log 165 TT 

64 ; 493018566815808 


-18974.368 


9 


, 1200329915 

-1925 vr log (220) tt 

^ ^ 9135984 


-33031.075 


10 


12155 , 325291539600149215255 

TT log (286) TT 

4 ^ 1172732412725203616 


-54866.421 


11 


, 31458443588344487293819 

-4641 TT log (364 ;:7r77:7:7:;r^r^^7r^;;;7r^:7:7^ 

^ ^ ' 60436675052957701680 


-87616.538 


12 


440895 25537984326378849719971131 
64 ^ ' 28270687046875000000000 ^ 


-135295.739 


13 


, 1779685691911133495 
-9975. log (560)- ,202109806542848 


-202952.031 


14 


56525 , 36234350694889865223938313068785 

TT log 680 TT 

4 ; 15613637127259094259005915136 


-296836.555 


15 


, 130243656594168370141034405 

— 19635 TT log (816) TT 

^ 37115886521993021558784 


-424587.139 



mentioned in Sec. [H they are needed together with those of the integer case in order to 
evaluate the information entropy of spherical and hyperspherical harmonics. Finally, it 
would also be desirable to extend the method introduced in this paper to other families 
of orthogonal polynomials having trigonometric representations, a line of research that 
is currently being developed. 
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Table 2. Exact and numerical values of the entropy E{C^^) for 1 < n < 15. 


n 


h;-v"q rti- Trti 1 n o 
JZjA.cLLL VcLiUc 


M 1 1 TVl OV1 /^Q 1 TTQ 1 1 1 O 

lilliciiCcLi VcLiUc 


1 


525 , 945 

TT log ( 5 I H TT 

128 s V ; -r ^Q24 


-17.839 


2 


2475 , , , 27685925 

— TT log (15) + ^^^^ TT 

128 ; ' 5225472 


-147.857 


3 


17325 , 61634724075 

TT log (35) + TT 

256 / i 3373232128 


-698.499 


4 


25025 , , 5573831525 
128 ' 115605504 


-2457.981 


5 


63063 , 338107973281463 

TT log (126) + TT 

128 ''^ ' 3173748645888 


-7150.909 


6 


143325 , 20887195 

77^7^ log (210) + TT 

128 ' 101376 


-18162.369 


7 


75075 , 1408430247274269205 

^^t:— TT log (330) H „^ . . tt 

32 ' 3944148534526464 


-41620.201 


8 


294525 , 806559968327725 

TT log (495) H TT 

64 / ' 1438588584576 


-87940.792 


9 


546975 , 29915266041851863399425 

TT log (715) H TT 

64 y . 37527437207206515712 


-173958.634 


10 


969969 , 97664804776687286561309 

TT log (1001) H TT 

64 ''^ ' 96698680084732322688 


-325775.232 


11 


6613425 , 230209361727271224010045 

TT log (1365) H TT 

256 ^ 212679240849405517824 


-582478.486 


12 


2723175 , , , 10188450005911283085 
64 log (1820)+ 12635587626401792 ^ 


-1000899.539 


13 


4352425 , 39663465263970548089202600252605 

TT log (2380) TT 

64 / 249818194036145508144094642176 


-1661590.212 


14 


6774075 , 717231543067734581588054629334401175 

IT log (3060) TT 

64 ^ 306761704661739640893688309874688 


-2676220.464 


15 


5148297 , 535111116210266542852402527915814650511 

TT log (3876) TT 

32 ^ 82233794352493419438828330115762176 


-4196611.889 
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